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SUMMARY 


A  solution  of  the  hydrodynamic  equations  is  presented  which  describes 
a  spherically  symmetric  shock  wave  in  a  polytropic  gas  resulting  from 
the  instantaneous  release  of  an  arbitrary  amount  of  energy  in  the  cen¬ 
ter.  The  external  force  of  gravity  as  well  as  some  variable  density 
and  pressure  distribution  of  the  undisturbed  medium  are  taken  into 
account.  It  is  shown  that  the  head  of  the  outwards  moving  shockwave 
is  located  at  a  distance  R(t)  =  0  •  from  the  center.  The  constant 
a  is  determined  from  the  adiabatic  exponent  of  the  gas.  the  density 
and  pressure  coefficients  of  the  ambient  medium  and  the  amount  of 
energy  supplied  by  the  explosion. 
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INTRODUCTION 


Sir  G.  I.  Taylor  [l "!  found  a  solution  of  the  hydrodynamic  equations 
which  represents  a  spherically  symmetric  strong -shock  wave  in  a 
polytropic  gas  resulting  from  a  very  intense  explosion  in  the  center. 

He  neglected  the  gravitational  effect  and  assumed  further  that  den¬ 
sity  and  pressure  of  the  ambient  medium  are  constant  and  zero, 
respectively.  The  latter  assumption  restricts  the  solution  automatic¬ 
ally  to  the  strong -shock  regime. 

Other  strong-shock  solutions  of  Taylor's  type  have  been  found  with 
regard  to  the  geometry  of  the  problem  [2',  and  with  regard  to  the 
inclusion  of  an  additional  heat  flux  term  in  the  hydrodynamic  equations 
[3"l.  The  spherically  symmetric  blast  wave  problem  including  the 
gravitational  effect  and  some  variable  density  distribution  of  the 
ambient  medium  has  been  solved  for  a  strong  shock  in  _4T. 

The  purpose  of  the  present  investigation  is  to  provide  a  spherically 
symmetric  solution  of  Taylor's  type  for  an  arbitrary  (strong  or  weak) 
shock  under  consideration  of  the  gravitational  effect  and  some  density 
and  pressure  distribution  of  the  ambient  gas  as  functions  of  the  dis¬ 
tance  from  the  center.  This  solution  may  be  applied,  for  example, 
to  the  description  of  the  hydrodynamic  phenomenae  resulting  from 
the  birth  of  a  nova. 
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ANALYTIC  FORMULATIONS  AND  ASSUMPTIONS 


The  problem  is  to  determine  radial  particle  velocity  u,  pressure  p 
and  density  0  within  an  outwards  moving  spherically  symmetric 
shock  wave  as  functions  of  radial  distance  r  from  the  center  and 
time  t.  In  the  case  of  a  polytropic  gas  with  the  adiabatic  constant  Y, 
these  quantities  have  to  satisfy  the  hydrodynamic  equations 


Sp  ,  1  9(r  ou) 

"  0 

c  t  i’  dr 

du  du  1  dp  _  d$ 

ol  dr  D  dr  dr 


a(ps  )  9(po  Y ) 

ht  Sr 


-  0 


(1) 


where  4  (r)  is  the  gravitational  potential  per  unit  mass.  A  further 
unknown  in  the  problem  is  the  position  r  -  R(t)  of  the  shock  front 
at  time  t  after  the  explosion  occurred. 


The  solutions  u,  p  and  0  of  the  hydrodynamic  equations  (1)  are  sub¬ 
ject  to  the  Rankine  -Hugoniot  transition  conditions  at  the  head  of  the 
shockwave.  L.et  ux ,  pa  and  p,  denote  the  values  of  u,  p  and  p 
immediately  behind  the  shock  front,  and  uD,  p0  and  P0  those 
immediately  ahead  of  the  shock  front.  The  Rankine -Hugoniot  rela¬ 
tions  are 

Pi  (R  -  u: )  -  0o  (R  -  u0)  ~  0 

Pi u j  (R  -  uj  -  c0uo.  (R  -  n0)  -  p1  -  p0  (2) 


3'-  (IUl3+  (R  •  u>}  ' 


-c  (t  u=s 


Y-l 


)  (R 


=  PlUl  -  Po  uo 

Here,  R  stands  for  the  velocity  of  the  shock  front. 
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The  second  restriction  on  u,  p,  and  o  is  that  the  amount  of  energy 
supplied  by  the  (instantaneous)  explosion  is  to  be  constant  with  res¬ 
pect  to  time  after  the  burst.  If  uo*(r),  po*(r)  and  Po  *(r)  denote 
the  particle  velocity,  pressure  and  density  of  the  gas  in  the  (station¬ 
ary)  case  prior  to  the  explosion,  this  amount  of  energy  is  given  by 

R(t) 

p  ___ 

E  =  4"T  j  I  y  (pu2  -  P^Uq*2)  -  (p  -  Qo*)  + 

^  ^  /Q\ 


—  (p  -  Pc*)  j  r2dr 

The  second  requirement  can  be  expressed  by  the  equation 


0,  t  >  to  =  0  (4) 

to  denotes  the  instant  when  the  explosion  occurs. 

We  assume  the  functions  uc*(r),  p0*(r),  p0*(r)  and  the  gravitational 
potential  ?(r)  to  be  of  the  form 

u*(r) 

P*(r) 

o*(r) 

Hr) 


=  0 


u 

■~3~ 

r 


3X_ 

h 


(5) 


dE  _ 
dt 

where  t  = 


Here  X  and  U  are  positive  constants.  Note  that  these  functions  do 
solve  the  hydrodynamic  equations  in  the  stationary  (or  more  pre¬ 
cisely:  in  the  quiescent)  case. 

For  r  =  R,  the  quantities  ue*(R),  po:|:(R)  and  0o*(R)  are  identical 
with  uc ,  p0  and  cu,  respectively.  The  Rankine -Hugoniot  relations 
(2),  therefore,  can  be  written  as 
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(Y+ 1)  R 


I  •  p  A. 

R  -  v  — 

L  u. 


R  j 


P- 


(Y+  1  )R  L 


2*a-  <Y-^'  k] 


(6) 


u(Y+j_) 

Ra 


r  _  R 

L  {Y- 1  )Ra  +  2v-l.  - 
u  R 


The  energy  supplied  by  the  explosion  is  then  given  by 
R(t) 

i*  r  i  3  3X  i  ,  p ,  i  , 

1  l7Pu  +  77“  ‘  “  (P  — s)  +  TpT  (P  ~ 

j  L  2  p  r  Y- 1 

0 


E  =  4TT 


r2  dr 


(7) 
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PROGRESSING  SHOCK  WAVE 


Following  a  classical  procedure,  we  try  to  solve  the  present  problem 
by  assuming  the  form  of  "progressing  waves"  for  the  solutions  in  order 
to  reduce  the  hydrodynamic  equations  to  a  system  of  ordinary  differ¬ 
ential  equations.  These  are  solutions  conveniently  written  as 


u(r ,  t)  =  a(r )  •  U(x) 

p(r,  t)  =  b(r )  '  P(x) 

o(r ,  t)  =  c (r )  ■  0(x) 

where  x  is  the  combination 
h(r) 

x(r,t)  = 


(8) 


(9) 


Inserting  equations  (8)  and  (9)  into  the  hydrodynamic  equations  (1), 
we  obtain  by  straightforward  computation  the  equations 


2rflhr 

4  W1  +  — 


r  ?ar  ~i  T 

|_(ac)'  +  ~  J  un+  ah'x  (Uf)'j  =  0 


a 

h 


{ 


-x3U '  +  a'hUs+  ah'xUU 1  + 


b'h 

ac 


P  bhd  P'|  _  3X_ 

Q  +  ac  X  111  "  '  U 


1 

(10) 


bc'Y  r  s 

—  Vx 


(PQ"Y),J-  U|  ah(^  -  Y^-)(PQ‘Y)  +  ah'x(Pa'Y)'jj-=  0 


Here  a  prime  denotes  the  derivative  with  respect  to  the  argument  of 
the  function.  Solutions  of  equations  (10)  can  be  obtained  by  separation 
of  the  variables  r  and  x.  The  separability  conditions  are 
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S 

n 


ah'  =  const, 
a'h  -  const, 
a  const. 


n  I  .  ‘.cil 

—  I  (ac) '  +  - 

c  _  r  -j 


const. 


b'h 

ac 

bh' 


(ID 


const. 


=  const. 


ah  |  -  Y  —  .  =  const. 

Lb  c  J 

The  first  three  equations  yield 


a(r) 


dr 


_-/a 


,  .  ,  I'*'?'* 

KM  --  .-j 


while  b  and  c  are  obtained  from  the  remaining  equations  as 
b(r)  -  3rJ)"1/s 
c(r )  =  Sr^ 

CL,  0,  6,o  and  Ware  constants. 

The  position  of  the  shock  can  be  determined  by  setting 

h(r)  /'r'N3^3  •  —  =  1  or  r  =  R(t)  =  Ot^3 

x  =  -  =  i  —  )  t 

t  Vcy 


( 1 2al 


(12b) 


(13) 
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Without  limitation  of  generality  we  can  choose 

a  =  1 

S  =  X  (14) 

5  -  u 


The  energy  supplied  by  the  explosion  is  then  given  by 


R(t) 

E  -  4r  !  |  -  u  r"‘ly<!  Q(x)U2(x)  +  3^  (r  Q(x)  - r”3  ) 

°o  L2 


+ 


X 

vTf 


sJU-  l/2 

(r 


P(x) 


r2  dr 


(15) 


In  order  that  the  energy  be  independent  of  time,  we  set 

^  =  .5 
2 


(16) 


Introducing  x  as  the  integration  variable  for  arbitrary  but  fixed  time  t, 
one  obtains 


E  = 


8-  p  r  1  ^  ? 

T  J  LimQ(x)u  (x)+  3x(n(x)-  1)  + 

"0 


X 

v  _  1 


(P(x)  -  l)j- 


dx 


[17) 


The  solution  for  the  shock  wave  is  represented  by  the  functions 


u(r,t)  = 

r‘l/s 

U(x) 

p(r,t)  = 

Xr’3 

P(x) 

C(r,t)  = 

Ur"2 

Q(x) 

X  = 

(■ -f 3 

'o' 

1 

t 

where  U(x),  P(x)  and  Q(x)  have  to  be  determined  from  the  system  of 
first  order  ordinary  differential  equations 
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(?.xs- 3£x  TT)  0'  +  £  (U-3XTJ*  )  fi  =  0 

'LM(2x2-  3exU)  U‘+£(^U3-6X)  jn  +  3\e(2P-xP')  =  0  (19) 


(2x3- 3€xU)(Pfi‘Y)'  +  2£  (2Y-3)  U  PH'7  =  0 


with 


e  =  (20) 

The  initial  conditions  for  U(x),  P(x)  and  fi(x)  are  obtained  from  the 
Rankine -Hugoniot  relations  at  the  shock  front  x  ~  1.  First  one  der¬ 
ives  from  equation  (13)  the  velocity  for  the  shock  front 

R  =  |ot'V3  =  JL  r'i/2  (21) 


Equations  (6)  become 


ui 


f  4u  -  9vXes-t  -l/3 

L3ia{v+l)e 


f8p  -  9(v-l)Xe3  ~ 
Pl  =  L  9{Y+l)e2  J 


(22) 


Pi 


f  2p2(y-tl) 
L2d(v_l)+9Y\€2 j 


On  the  other  hand,  the  same  values  for  ux  ,  p1  and  px  are  obtained 
from  equations  (18)  if  we  set  r  =  R  and  x  =  1,  The  desired  initial 
conditions  are  therefore 


U(l) 


4M-9YXe3 

"3)j(Y+l)e 


(23) 
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P(l) 


8U  -  9(Y-l)>-e2 
9(Y+1)U 


rf(  i) 


2m(Y+1) 
2i-i(v-1)+9v>  ^ 


(23) 


The  integration  of  the  differential  equations  (19)  can  be  reduced  to 
the  problem  of  solving  a  second  order  ordinary  differential  equation 
for  U(x)  and  determining  P(x)  and  0(x)  by  quadratures.  With  the 
abbreviations 


£  ( 3xU  '  -  U  )  _  $ 

2x2-  3£xU 

2£(2Y-3)U  _ 

2x2-3exU 

Try—  I  2x“  U '  -  e(3xU 1  -  U)  U  i  -2  ::  A 

J  At  L  J 

it  follows  from  the  first  and  third  equation  in  (19)  that 

^x 

J  $dx 

fi(x)  -  il(l)  e1 

x 

J  (V+Y$)  dx 

P(x)  ”  P(l)  e1 


(24) 


(25) 


These  explicit  expressions  for  ^(x)  and  P(x)  lead  to 


xP'-2P 


r': X '+4./V  1 1  a>  Cx(*+y$)-2T  ;  . 

,1.  ei’*  X.  "  x(Y+v$).2  J  x 


n  Q( 

The  second  equation  in  (19)  reads 


A  i 

iP'  -2P  ?  .  1J  A”  dx 

=  a  =  -'(I)  e 


(27) 


Q 
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Equations  (26)  and  iZi)  furnish 


Pjl) 

^d)A(i) 


d'i'(l)av§(l)-2]e1'’ 


|  Y+(v 


-1)$+ 


[  x(v-‘-v$)-2ll 
x(7+y$)  -2 


-7 —  P  dx 

i\  J 

(28) 


This  equation  holds  identically  in  x,  if  and  only  if  the  equations 


Y +  (y. 1) $  + 


[x(Y+ Y$) -2] ' 

x(Y+ y?)  -2 


A* 


0 


P(l) 

Q(1)A(1) 


[Y(l)+  v§(l)  .  2] 


1 


(29) 


are  satisfied.  Here  the  first  equation  represents  the  second  order 
ordinary  differential  equation  for  U(x),  while  the  second  equation 
furnishes  the  required  initial  value  U'(l)  on  substitution  of  (24). 

For  given  physical  constants  Y,  X  and  M,  the  solutions  of  equations 
(25)  and  (29)  will  be  functions  of  x  and  of  the  parameter  £.  If  Y  =  3/2, 
the  flow  is  isentropic. 

It  is  meaningful  to  ask  under  what  conditions  the  pressure  p  and  the 
density  P  remai n  positive.  Equations  (18)  require  that  P(x)  and 
(x)  be  positive.  According  to  equations  (25),  this  is  the  case  if 
the  integrals  over  ?  and  ¥  exist  and  both  constants  Q(  1 )  and  P(l) 
are  positive.  It  follows  then  from  the  second  of  equations  (23)  that 
the  parameter  S  has  to  be  restricted  to  the  interval 


2 

0  <  e  <  - 


2u 


3  ( Y  - 1 )  A 


(30) 


Equation  (17),  after  inserting  the  solutions  U(x,  e),  P(x,  e)  and 
(x,  e)  and  integrating,  gives  the  relation  between  the  energy  E 
supplied  by  the  explosion  and  the  parameter  S.  If  £  varies  between 
the  limits  given  in  (30),  E  may  range  between  zero  and  infinity. 
Solving  equation  (17)  for  £  (or  tJ),  one  arrives  finally  at  the  result 
that  the  position  of  the  shock  front  at  time  t  after  the  explosion  is 
determined  by  the  distance 

r  =  R(t)  -  o(v,  ,\,u,  E)  •  t3  (31) 
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It  is  interesting  to  compare  this  result  with  that  obtained  by  Taylor 
C  1 1  in  the  case  of  a  strong  shock,  where  density  and  pressure  out¬ 
side  the  wave  are  assumed  to  be  constant  and  zero,  respectively, 
and  gravity  is  neglected.  His  result  states  that  the  shock  front  can 
be  localized  by  the  distance  r  -  R(t)  ~  CJ*  •  ,  where  the  constant 

O  *  is  determined  in  a  similar  manner  as  3  from  the  adiabatic 
constant  of  the  gas,  the  outside  (constant.)  density  and  the  energy 
supplied  by  the  explosion. 


Comparison  of  the  two  results  shows  that  the  expansion  velocity  of 
the  wave  in  the  present  case  is  greater  than  that  of  Taylor's  wave, 
at  least  from  a  certain  time  on,  depending  on  the  ratio  of  the  constants 
o  and  cpN.  This  fact,  of  course,  follows  from  the  different  assumptions 
made;  it  might  be  accounted  for  by  the  assumed  density  and  pressure 
distribution  only,  since  the  gravity  might  be  expected  to  cause  the 
opposite  effect.  Nevertheless,  it  is  remarkable  that  this  is  true  also 
if  the  shock  is  weak. 
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